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In  order  to  offer  a  guideline  to  search  for  new  oxyapatites  exhibiting  high  ionic  conductivities  (or)  for  possi¬ 
ble  application  in  solid  oxide  fuel  cells  operating  at  intermediate  temperatures  (IT-SOFCs),  a  mathematical 
model  is  developed,  i.e.,  a  linear  relation  of  a  and  the  relative  Coulomb  energy,  to  predict  the  a  of  oxyap¬ 
atites.  The  inputs  to  the  model  are  two  readily  available  elemental  properties,  namely  the  ionic  radius  and 
the  electronegativity  of  the  constituent  elements.  The  model  not  only  predicts  the  ionic  conductivities 
of  45  oxyapatites  but  also  rationalizes  the  observed  trends  reported  in  the  literatures.  The  effects  of  the 
two  elemental  properties  on  ionic  conductivity  are  also  examine  and  predict  a  wide  range  of  new  oxygen 
stoichiometric  oxyapatites  with  ionic  conductivities  potentially  as  high  as  10~2-10-1  Scm-1  at  500  °C  is 
predicted.  The  investigations  suggests  an  optimization  strategy  to  search  for  promising  oxyapatites,  i.e., 
applying  dopants  with  large  ionic  radii  and  low  electronegativities. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Materials  with  high  ionic  conductivities  have  attracted 
widespread  interest  due  to  their  possible  application  in  many  areas 
such  as  solid  oxide  fuel  cells  (SOFCs),  oxygen  sensors  and  separa¬ 
tion  membranes  [1-4].  Yttria  stabilized  zirconia  (YSZ)  is  the  most 
favoured  SOFC  electrolyte  due  to  its  high  oxygen  ionic  conductiv¬ 
ity  and  structural  stability  at  elevated  temperatures  (900-1000  °C). 
On  reducing  the  operating  temperature  of  a  SOFC  to  intermediate 
values  of  600-800  °C  (IT-SOFCs),  however,  the  ionic  conductivity 
of  YSZ  decreases  significantly  [2].  Thus,  to  overcome  the  disadvan¬ 
tages  arising  from  a  reduced  operating  temperature,  great  efforts 
have  been  made  to  search  for  other  oxygen  ion  conductors  with 
high  ionic  conductivity  at  lower  temperatures  (600-800  °C)  [1,5]. 
Among  them,  an  apatite-type  oxygen  ion  conductor  based  on 
lanthanum  silicate,  La9.33Si6026,  shows  significant  promise  as  an 
electrolyte  for  of  IT-SOFCs  [1-19].  It  has  been  reported  [4]  that  the 
conductivities  of  La9.33Si6026  are  1.1  x  10-4  and  2  x  10-3  Son-1  at 
500  and  800  °C,  respectively. 
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The  structure  of  a  typical  oxyapatite  consists  of  isolated  Si04 
tetrahedral  units  and  rare  earth  cations  located  in  9  coordinate 
(Lal/La2)  and  7  coordinate  (La3)  channel  sites  [20].  The  remain¬ 
ing  oxygen  ions  occupy  one-dimensional  channels  running  through 
the  structure,  which  are  considered  as  vital  for  high  oxygen  ion  con¬ 
duction  in  such  materials  [21].  Both  experimental  and  theoretical 
evidence  indicate  that  oxyapatites  conduct  ions  via  an  intersti¬ 
tial  mechanism  [1,7,22],  and  local  cooperative  relaxation  of  Si04 
tetrahedra  is  an  important  factor  for  facile  oxygen  ion  conduc¬ 
tion.  In  other  words,  Si  and  the  nearest  oxygen  ions  generally 
displace  away  from  the  channel  and  towards  the  Lal/La2  sites,  as 
suggested  by  the  observed  changes  in  ionic  conductivity  of  Mg- 
doped  La9.33(Si04)602  [23].  A  large  enhancement  in  conductivity 
was  observed  when  Mg  was  doped  on  the  Si  sites,  whereas  the 
conductivity  decreased  when  Mg  was  doped  on  La  sites  [23].  This 
suggests  that  Mg  leads  to  a  local  expansion  of  the  channels  when  it 
is  on  Si  sites,  while  it  encroaches  the  oxygen  ion  channels  when  it  is 
doped  on  La  sites.  Further  studies  revealed  that  non-stoichiometry 
and  partial  occupancy  on  the  Lal/La2  sites  have  the  same  effect 
on  the  conductivity  as  the  incorporation  of  excess  oxygen  into  the 
system  [8].  Also,  it  was  reported  that  a  wide  range  of  cation  sub¬ 
stitutions  in  an  oxyapatite  are  possible  due  to  the  flexibility  of  the 
apatite  structure  in  accommodating  a  range  of  ion  sizes,  and  this 
enhances  the  success  in  searching  for  new  oxyapatites  of  high  ionic 
conductivity. 


Y.  Zeng  et  al  /  Journal  of  Power  Sources  196(2011)  4524-4532 


4525 


Inspired  by  the  above  findings  from  experimental  and  theoreti¬ 
cal  investigations,  a  mathematical  model  is  proposed  in  this  study 
to  predict  the  ionic  conductivities  of  oxypaptites  in  an  attempt  to 
provide  a  guideline  to  search  for,  or  synthesize,  new  oxyapatites. 
As  an  initial  attempt,  only  oxygen  stoichiometric  oxyapatites  are 
considered.  In  the  model,  it  is  postulated  that  the  expansion  of  c- 
axis  channels  plays  a  central  role  in  the  enhancement  of  oxygen 
ion  conduction  and  the  Coulomb  interaction  between  two  charged 
units  of  Si04  tetrahedra  and  a  Lai  /La2  site  dominates  the  expan¬ 
sion  or  contraction  of  c-axis  channels.  The  interaction  energy,  i.e., 
the  relative  Coulomb  energy,  is  formulated  by  a  set  of  variables 
that  are  properly  defined  in  terms  of  two  fundamental  properties 
of  the  constituent  elements  (i.e.,  electronegativity  and  ionic  radius) 
and  composition.  According  by,  a  linear  relation  between  the  rel¬ 
ative  Coulomb  energies  and  the  logarithmic  ionic  conductivities  is 
established.  This  linear  relation  is  examined  for  45  oxygen  stoi¬ 
chiometric  oxyapatite  compositions  reported  in  the  literature,  and 
a  good  correlation  with  the  coefficient  of  determination  R 2  =  0.81  is 
achieved.  Since  the  inputs  to  the  model,  namely,  the  fundamental 
properties  of  elements,  are  readily  available,  the  model  enables  the 
evaluation  or  prediction  of  the  conductivities  for  a  large  number  of 
oxyapatites  with  much  less  computation  effort.  The  effects  of  ionic 
radius  and  electronegativity  on  ionic  conductivity  are  discussed  on 
the  basis  of  the  model.  An  examination  is  also  made  of  oxyapatites 
that  are  virtually  generated  by  numerically  substituting  a  range  of 
elements  into  the  Si  and  Lai  / La2  sites,  respectively.  New  oxyap¬ 
atites  with  high  conductivities  from  10-2  to  10-1  Scm-1  at  500  °C 
are  predicted. 


2.  Model  development 

Oxygen  stoichiometric  apatite  oxides  have  the  general  formula 
Aio-x(B04)602  where  A  and  B  are  a  range  of  cations,  e.g.,  A  is  a 
rare-earth  or  alkaline-earth  element,  and  B  is  Si  or  Ge,  etc.  The 
structure  of  the  typical  oxyapatite  La9.33Si6026  is  hexagonal  in  sym¬ 
metry  where  the  B  cations  (Si)  take  tetrahedral  coordination  with 
the  structural  oxygen  (Si04),  whereas  the  A  cations  (La)  locate  in 
seven-  (La3)  and  nine-coordinate  (Lal/La2)  cavity  sites.  The  nine- 
coordinate  position  lies  in  the  centre  of  the  smaller  of  the  two 
structural  channels,  and  the  seven-coordinate  position  lies  around 
the  periphery  of  the  larger  channel  in  the  centre  of  which  reside 
oxygen  ions.  The  extra  oxygen  ions  occupy  the  larger  channels  run¬ 
ning  along  the  c-axis,  which  are  responsible  for  the  high  oxygen  ion 
conductivity  [14-16]. 

As  stated  in  [24],  lanthanum  atoms  usually  fully  occupy  La3 
sites,  while  cationic  vacancies  are  randomly  distributed  on  Lai  /La2 
sites  where  rare-earth  and  alkaline-earth  elements  can  be  doped. 
Therefore,  the  interaction  between  the  Si04  tetrahedra  and  Lai  /La2 
site  is  postulated  to  be  dominant  for  the  expansion  or  contraction 
of  the  c-axis  channel.  In  this  study,  the  interaction  is  represented 
by  a  relative  Coulomb  energy  between  the  two  charged  sites  of  A 
(Lal/La2)  and  B  (Si04),  as  illustrated  in  Fig.  1.  On  the  other  hand, 
the  Coulomb  interaction  between  the  Si04  tetrahedra  and  a  La3  site 
is  secondary  since,  in  most  cases,  the  La3  site  remains  intact  when 
the  chemistry  of  an  oxyapatite  is  changed.  In  another  words,  the 
Si04  tetrahedra/La3  interaction  forms  a  common  reference. 

The  A  and  B  sites  are  taken  as  two  circles  as  shown  in  Fig.  2. 
The  relative  Coulomb  energy  between  the  two  sites  is  formulated 
by  two-dimensional  integration  of  the  charge  densities  of  the  two 
sites.  The  derivation  starts  from  the  classic  Coulomb  energy,  as 
defined  by  Eq.  (1)  where  qA  and  qg  are  charges  of  A  and  B  sites, 
respectively,  and  r  and  ke  denote  the  distance  between  the  centres 
of  the  two  sites  and  the  Coulomb  constant,  respectively. 


Fig.  1.  Two  charged  units  in  structure  of  La9.33Si6  026  (viewed  from  c-axis): 
(A)— Lal/La2;  (B)— tetrahedral  (Si04). 
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The  Coulomb  energy  is  evaluated  by  assuming  that  the  charge 
densities  of  the  two  units  are  uniformly  distributed  along  the 
peripheries  of  A  and  B,  respectively.  The  charges  of  two  infinitesi¬ 
mals  of  circles  A  and  B,  i.e.,  AA  and  AB,  are  evaluated  as  Eqs.  (2)-(3) 
where  pA  and  pB  are  the  charge  densities  of  the  peripheries  of  A  and 
B,  and  RA  and  R B  are  the  ionic  radii  of  A  and  B,  respectively.  It  is  noted 
that  RB  is  the  sum  of  the  radii  of  Si4+  and  O2-. 


A  qA  =  pARAd6A 

(2) 

AqB  =  pBRBdOB 

(3) 

The  distance  between  the  two  infinitesimals  on  circles  A  and  B, 
denoted  as  r,  is  evaluated  as: 

r2  =  c(l  +  a  sin  0B  +  b  cos  Ob) 

(4) 

where 

—2 RaRb  sin  0A 

°  ~  Rj  +  R2B  +  d2  -  2 dRA  cos  0A 

(5) 

^  _  -2Rs(d  -  RA  cos  Pa) 

~  R2A  +  R2B  +  d2  -2dRA  cos  0A 

(6) 

c  =  RA+R-l+d2  -  2dRA  cos  0A 

(7) 

The  Coulomb  energy  between  circle  A  and  circle  B  can  be  eval¬ 
uated  by  taking  the  double  integrations  as  Eq.  (8). 


Ea-b  =  PaPb^aRb 


l*2jt  i*2jt  ^ 

Jo  Jo  T 


-dOBdOA 


Substituting  r  (defined  in  Eq.  (4))  into  Eq.  (8)  yields: 


Ea-b  =  PaPbRaRb  /  -p 


n2n  1  p2n-\-a 

l 


1 


y/\  +  sin  0 


dO 


(8) 


(8a) 


where  a  =  arctan(a/b)  and  0  =  a  +  0B.  Let  y  =  jJ7r+“  1  /yf  1  +  sinOdO 
and  note  that  y  turns  out  to  be  a  constant.  Substitution  of  c  and  y 
into  Eq.  (8a)  gives: 


Ea~b  = 


YPaPbRaRb 


n2l T 

JO 


^/R2+R2+d2  JQ  0^ 

cos  0A 


dOA 


where 


—2dRfi 

+  R2  +  d2 


(8b) 


(9) 


Ea-b  = 


Let  Z  =  ff  \l \J\  +  §  cos  0A6QA  and  then  Eq.  (8b)  is  reduced: 
Y^PaPbrarb 


y/Rl+RZ+V 


(8c) 


In  this  work,  the  charge  density  p  and  the  circle  radius  are  asso¬ 
ciated  with  the  electronegativities  and  ionic  radii  of  the  constituent 
elements,  respectively.  Both  electronegativity  and  ionic  radius  are 
fundamental  elemental  properties.  The  resultant  Coulomb  energy 
Ea_b  is  thus  a  relative  quantity.  To  establish  the  relationship 
between  the  relative  Coulomb  energy  and  the  ionic  conductivity 
with  acceptable  accuracy,  the  charge  densities  and  ionic  radii  in  Eq. 
(8c)  are  re-defined  with  the  properties  of  La  and  Si  as  the  reference 
in  the  following  section. 

The  general  formula  of  an  oxygen  stoichiometric  oxyapatite 
with  dopants  is  given  as  Al alA2a2Bl blB2b2B3b3026  where  the  A  and 
B  sites  are  occupied  by  up  to  two  and  three  types  of  cation,  respec¬ 
tively.  In  particular,  Al  is  occupied  by  trivalent  cations  only  and  A2 
can  be  occupied  by  divalent  or  trivalent  cations;  B\  is  occupied  by 
tetravalent  cations  only,  while  B 2  and  B3  can  be  occupied  by  either 
divalent  or  trivalent  cations. 


Table  1 

Elemental  properties  of  ions. 


Element 

Charge 

Ionic  radius  (A) 

Electro-negativity 

Al 

3+ 

0.54 

1.61 

B 

3+ 

0.23 

2.04 

Ba 

2+ 

1.35 

0.89 

Ca 

2+ 

1 

1.00 

Co 

3+ 

0.55 

1.88 

Eu 

3+ 

0.95 

1.12 

Fe 

3+ 

0.55 

1.83 

Ga 

3+ 

0.62 

1.81 

Gd 

3+ 

0.94 

1.20 

Ge 

4+ 

0.53 

2.01 

La 

3+ 

1.03 

1.10 

Mg 

2+ 

0.72 

1.31 

Nd 

3+ 

0.98 

1.14 

0 

2- 

1.4 

3.44 

Pr 

3+ 

0.99 

1.13 

Si 

3+ 

0.4 

1.90 

Sm 

3+ 

0.96 

1.17 

Sr 

2+ 

1.18 

0.95 

Zn 

2+ 

0.74 

1.65 

The  subscripts  al,  a2,  hi,  b2  and  b3  designate  the  composi¬ 
tions  of  the  associated  elements,  and  the  total  number  of  cations 
is  summed  up  as  N= al  +  a2  +  b\  +  b2  +  b3.  The  average  charges  of  A 
and  B  sites  are  calculated  as  Eqs.  (10)  and  (11),  respectively. 


n™g  =  xnA  i  +(1  -x)nA2 

(10) 

nBvg  =y\nB\  +y2«B2  +  (i  -yi  -y2)nB3 

(11) 

where  nA1 ,  nA2,  nB  i ,  nB2  and  nB3  are  the  valences  of  cations  Al  ,A2,  B\ , 
B2  and  B3,  respectively;  x,yi,  andy2  are  normalized  compositions, 
i.e.,  x  =  al/(al  +a2),  y1  =  bl/(M  +b2  +  b3),  y2  =  b2\(b\  +b2  +  b3).  The 
superscript  ‘avg’  denotes  the  average  values. 

The  average  ionic  radius  of  the  A  site  is  evaluated  in  terms  of  the 
radii  of  Al  and  A2.  Since  Al  can  only  be  occupied  by  trivalent  cations, 
its  radius  remains  unchanged  and  is  denoted  by  the  superscript 
‘org\  On  the  other  hand,  A2  can  be  occupied  by  divalent  or  trivalent 
cations,  its  radius  is  thus  modified  as  Eq.  (12)  and  denoted  by  the 
superscript  lmod\  Specifically,  the  A2  original  radius  is  multiplied 
by  a  factor  defined  as  the  valence  ratio  of  A2  to  La3+  (=3). 


nmod  _  ™A2  p org 
*712  -  3  nA2 


(12) 


The  average  radius  of  the  A  site  is  defined  as: 

Rafg  =  xR°Jg  +  ( 1  -  x)R™od  (13) 


The  ionic  radii  of  B2  and  B3  are  modified  in  terms  of  the  Si 
valence,  nsi  (=4),  and  the  coordination  numbers  of  Si4+  (=4)  and 
La3+  (=9). 


nmod 

KB2 


nB2 

nsi 


x 


9 

4 


xC 


(14) 


nmod 

KB3 


nB3 

"si 


9 

4 


(15) 


The  original  radius  of  B\  is  used  without  any  modification  as  B\ 
can  only  be  occupied  by  tetravalent  cations.  The  average  radius  of 
the  B  site  is  thus  obtained  in  the  similar  fashion  as  that  of  the  A  site, 
namely: 

R7s=yKf  +y2Rf2d  +  n  -yi  -y2)K£od  06) 


The  distance  between  the  centres  of  A  and  B  sites  can  be  repre¬ 
sented  as  the  sum  of  the  average  radii  of  the  A  site  the  and  B  site, 
multiplied  by  a  constant  X  (>0),  i.e., 


d  =  X(Ra;s  +  Rlvs  +  R°0rs) 


(17) 
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Table  2 

Reported  oxyapatite  ionic  conductivities  at  500  °C. 


No. 

Oxyapatite 

Conductivity,  0  (S  cm-1 ) 

References 

1 

Prg.33Si6026 

8.10  X  10-5 

Table  1  in  [2] 

2 

Lag.33Si6026 

1.10  X  10-4 

3 

Ndg.33Si6026 

1.00  X  10-4 

4 

Smg.33Si6026 

2.20  X  10-5 

5 

Gdg.33Sie026 

1.50  X  10-6 

6 

Las.67BaSi6026 

1.40  X  10-4 

7 

La8.33Bai.5Si6026 

6.60  x  10-5 

8 

Lag.67SrSi6026 

8.30  x  10-5 

9 

La8.33Sri.5Si6026 

1.30  x  10-5 

10 

La8.67CaSie026 

5.80  x  10-5 

11 

La8.33Cai.5Si6026 

3.40  x  10-5 

12 

Lag.5Si5.5Gao.5O26 

4.60  x  10“4 

13 

Lag.67SisGa026 

1.00  x  10“3 

14 

Lag.83Si4.5Ga1.5O26 

1.30  x  10-3 

16 

Lag.5Si5.5B0.5O26 

4.10  x  10“4 

17 

Lag.67SisB026 

3.30  x  10-4 

18 

Lag.83si4.5Bl.5026 

4.90  x  10-4 

19 

LagMgo.5Sie026 

2.10  x  10-5 

Table  2  in  [9] 

20 

Las.67MgSi6026 

3.60  x  10-6 

21 

Lag.33  Mgo.5  Sis.s  Mgo.5  O26 

1.60  x  10-4 

22 

Lag.5  Si5.75  Mgo.25  026 

1.80  x  10“3 

23 

Lag.67si5.5Mgo.5026 

3.00  x  10-3 

24 

Lag.5  Si5.75  Zno.25  O26 

2.10  x  10-3 

25 

Lag.33  Bao.5  Sis.5  Mgo.5  026 

1.80  x  10-3 

Table  3  in  [9] 

26 

Lag. 83  Sis  Mgo.5  Bo.5  O26 

1.90  x  10-3 

27 

Lag. 83  Sis  Mgo.5  Gao.5  026 

3.10  x  10-3 

28 

EUg.33Si6026 

1.30  x  10-5 

Table  4a  in  [9] 

29 

Prg.67SisGa026 

1.10  x  10-4 

Table  4b  in  [9] 

30 

Ndg.67SisGa026 

1.50  x  10-5 

31 

Smg.67Si5Ga026 

7.60  x  10-6 

32 

Eu9.67Si5Ga026 

7.20  x  10-6 

33 

Gdg.67SisGa026 

1.20  x  10“6 

34 

Lag.83si4.5Fel.5026 

1.26  x  10-4 

Fig.  10  in  [11] 

35 

La6.83  Pf3  Si4.s  Fei.5  O26 

1.00  x  10-4 

36 

La3.83pr6si4.5Fel.5026 

1.00  x  10-4 

37 

Lag.83si4.5All.5026 

1.00  x  10“3 

38 

Lag. 83  Si4.5AlFe0.5  O26 

2.51  x  10“4 

Fig.3  in  [10] 

39 

Lag.83Si4.5Alo.5Fe026 

1.26  x  10-4 

40 

Lag.67SisFe026 

1.03  x  10-4 

Fig.  2  in  [8] 

41 

Lag.83Si4.5CO1. 5  O26 

8.16  x  10-4 

Fig.  4  in  [8] 

42 

Lag.67  Sis  C0O26 

5.15  x  10-4 

43 

La8.5srl.5si5.5Alo.5026 

8.20  x  10-6 

Table  1  in  [12] 

44 

Lag.5si5.5Alo.5026 

1.70  x  10“4 

45 

Lag.33  Ge6026 

1.12  x  10-4 

Fig.  4  in  [13] 

Eq.  (9)  can  thus  be  rewritten  as: 

-2dRa;s 

^Cs)2+(Cs+Cs)2  +  d2 


(9a) 


Next,  the  electronegativities  of  the  A  and  B  the  sites  are  modified 
as  represented  by  Eqs.  (18)-(19),  where  E^8  (M=A1,  A2,  B\,  B2,  B3) 
denotes  the  original  electronegativity  of  the  corresponding  cation. 


EAmod  =  *E°XS)06  +  (1  -  *)EM20d)°'6 


(18) 


Erd  =  y^sW)0'6  +y2EB?(RB2od)0'6 

+  (1  -yi  -y2)£B3S(RSod)°'6  (19) 


Similarly,  the  modified  electronegativities  of  02~  and  La3+  are 
evaluated  as  Eqs.  (20)  and  (21),  respectively. 

E™d  =  E00rs(R°rsf6  (20) 


Fig.  3.  Linear  relationship  between  experimental  ionic  conductivities  (in  loga¬ 
rithmic  form)  and  calculated  nominal  Coulomb  energy  of  oxygen  stoichiometric 
oxyapatites  reported  in  literature. 


EL7d=C(Ru)06  (21 ) 

The  charge  densities  of  the  A  site  and  B  site,  i.e.,  p™od  and  p™od, 
are  defined  in  terms  of  the  modified  electronegativities  and  the 
average  charges,  as  shown  in  Eqs.  (22)  and  (23),  respectively.  In 
particular,  p™od  is  defined  as  a  standard  density  of  unity  minus 
a  deviation  from  the  electronegativity  of  La3+  when  the  A  site 
is  substituted  with  other  cations.  The  deviation  is  normalized  by 
the  modified  electronegativity  of  La3+  (£™d)  and  averaged  by  the 
charge  of  A  site  (n“vg).  p™od  is  also  defined  in  the  similar  way  where 
the  deviation  from  the  standard  density  of  unity  is  caused  by  the 
deviation  of  the  B  site  from  O2-.  It  is  noted  that  the  normalized 
deviation  is  multiplied  by  the  factor  of  (9/4)2  with  9  and  4  being 
the  coordination  numbers  of  Lai  /La2  and  Si4+,  respectively.  Since 
positive  Coulomb  energy  is  taken  into  account,  the  absolute  value 
of  p™od  is  used  as  shown  in  Eq.  (23). 


cmod  _  c-mod 
muod  _  i  _  nA  nLa 

c-mod  navg 

cLa  nA 


(22) 


nmod 

Pb 


-1  + 


c-mod  _  c-mod 
nB  nO 
cmod  „Eivg 

Eq  lift 


(23) 


The  re-defined  radii  in  Eqs.  (13)  and  (16)  and  charge  densities 
in  Eqs.  (22)-(23)  are  substituted  into  Eq.  (8c).  By  removing  the 
constant  y,  the  relative  Coulomb  energy  is  obtained  as  defined  by: 


E-mod 

CAB 


Zp™Ap™AR™g{R™s  +  R°rg) 

V(«T)2+(«?*+02+^b 


(8d) 


Fig.  4.  Comparison  of  ionic  conductivities  of  oxyapatites  NdxSiyGa(6_y)026  obtained 
from  model  prediction  (diamonds)  and  literature  (squares  from  [24]  and  asterisks 
from  [2,9]). 
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Fig.  6.  Calculated  ionic  conductivity  by  substituting  one  tetravalent  cation  at  Si  site. 
General  formula  of  oxyapatite  is  La9.33Si5M026  (M  =  M4+). 


universal  constant  X,  as  introduced  in  Eq.  (17),  is  given  as  1,  imply¬ 
ing  close  contact  of  two  charged  units. 

F^od  is  calculated  and  divided  by  the  total  number  of  cations  for 
each  sample  listed  in  Table  2,  and  plotted  against  the  logarithmic 
ionic  conductivity  (Fig.  3).  A  linear  relation  is  found  with  the  coeffi¬ 
cient  of  the  determination  ( R 2 )  of  0.81 ,  which  indicates  a  very  good 
linear  correlation. 

The  linear  function  is  used  to  predict  the  ionic  conductivity  from 
the  relative  Coulomb  energy.  The  predictive  capability  is  exam¬ 
ined  by  comparing  the  calculated  ionic  conductivities  with  the 
reported  results  [2,9,27]  of  five  oxyapatites  of  the  general  formula 
NdxSiyGa(6_y)026  (y  =  4-6),  as  shown  in  Fig.  4.  The  results  show 
that  the  ionic  conductivities  of  these  oxyapatities  decrease  as  more 
Ga3+  cations  are  doped  at  the  Si4+  site,  and  the  predicted  trend  is 
in  agreement  with  the  literature  data.  While  the  ionic  conductiv¬ 
ities  generated  from  the  model  decrease  monotonously,  the  data 
reported  in  the  literature  fluctuate  around  the  trend  line,  with  the 
exception  of  that  for  NdioSi4Ga2026. 

3.2.  Effects  of  doping 


Fig.  5.  Calculated  ionic  conductivity  (Scm-1)  by  doping  M  at  La  and  Si 
sites,  (a)  M  =  Mg2+,  Ca2+,  Sr2+,  and  Ba2+  (shaded  circles— La8.67MSi6  026, 
open  circles— La8M2Si6  026,  shaded  diamonds— LaioSi5M026,  open 

diamonds— La9.67Si5.5M0.5O26);  (b)  M  =  A13+,  Ga3+,  In3+  and  Tl3+ 

(shaded  circles— Las.33MSi6026,  open  circles— La7.33M2Si6  026,  shaded 
diamonds— La9.67Si5M026,  open  diamonds— La9.5Si5.5M0.5O26);  (c)  M  =  Bi3+, 

Nd3+,  Pr3+  and  Ce3+  (shaded  circles— La8.33MSi6026,  open  circles-La7.33M2Si6  026, 
shaded  diamonds  La9.67Si5M026,  open  diamonds— La9.5Si5.5M0.5O26). 


The  relative  Coulomb  energy  F™d  is  then  correlated  to  the  ionic 
conductivity  <7.  As  shown  in  the  following  section,  F™d  is  nearly  lin¬ 
ear  to  the  logarithmic  <7.  It  follows  from  the  model  that  the  larger 
the  relative  Coulomb  energy,  the  higher  the  ionic  conductivity.  Fur¬ 
thermore,  it  is  readily  available  to  predict  the  ionic  conductivity  of 
oxyapatite  electrolyte  materials  by  using  this  model. 

3.  Results  and  discussion 

3.1.  Model  validation 

As  shown  in  Eq.  (8d),  the  inputs  to  the  calculation  of  £™d  are 
the  fundamental  element  properties  of  ionic  radii,  electronegativ¬ 
ities,  and  valences  that  can  be  found  in  [25,26].  Table  1  lists  the 
fundamental  data  of  the  elements  used  in  this  study.  The  ionic  con¬ 
ductivities  of  45  reported  oxyapatites  at  500  °C  are  summarized  in 
Table  2  and  are  used  to  validate  the  developed  model.  The  unknown 


In  this  section,  the  effects  of  doping  cations  on  the  changes  of 
ionic  conductivity  are  investigated  based  on  the  typical  oxyapatite 
La9.33Si6026  by  using  the  developed  model.  In  particular,  ionic  con¬ 
ductivities  predicted  from  doping  cations  of  valence  2+,  3+  or  4+  on 
the  La  or  Si  sites  are  examined.  Since  the  present  study  focuses  on 
oxygen  stoichiometric  oxyapatites,  the  stoichiometry  of  the  cations 
is  adjusted  accordingly  to  maintain  the  oxygen  stoichiometry  at  26 
when  incorporating  dopants  into  the  structure.  In  order  to  conduct 
a  consistent  study  among  all  the  predictions,  the  predicted  ionic 
conductivity  a  (7.36  x  10-4  S  cm-1 )  instead  of  the  experimental  a 
(1.1  x  10_4Scm_1)  of  La9.33Si6026  is  taken  as  the  reference. 


Fig.  7.  Calculated  ionic  conductivity  of  oxyapatites  La9.67Si5M026  (M  =  3+)  and 
La9.33Si5M026  (M  =  4+)  by  doping  M  cations  of  4th  period  (Ni,  Co,  Fe,  V,  Ti)  and  5th 
period  (Pd,  Rh,  Ru,  Nb,  Zr). 
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3.2 A.  Doping  at  the  La  site 

The  effects  of  doping  cations  on  the  La  site  of  Lag.33Si6026  are 
studied  with  three  simulated  scenarios  that  correspond  to  doping 
three  groups  of  cations,  respectively,  namely:  (i)  Mg2+,  Ca2+,  Sr2+, 
and  Ba2+  (Group  IIA  in  the  periodic  table);  (ii)  Al3+,  Ga3+,  In3+,  and 
Tl3+  (Group  IIIA  in  the  periodic  table);  (iii)  Bi3+,  Nd3+,  Pr3+,  and  Ce3+. 

Doping  divalent  cations  of  Group  (i)  at  the  La  site  was  first  stud¬ 
ied.  Specifically,  the  incorporation  of  two  levels  of  dopants  M  (Mg2+, 
Ca2+,  Sr2+,  and  Ba2+)  was  investigated  according  to  the  formula 
La8.67MSi6026  and  La8M2Si6026.  It  is  noted  that  the  lower  charge 
on  M  (2+)  is  balanced  by  an  increased  La  content  so  that  the  oxy¬ 
gen  content  is  maintained  as  nominally  stoichiometric.  As  shown 
in  Fig.  5a  (shaded  and  open  circles),  all  the  calculated  ionic  conduc¬ 
tivities  (or)  are  smaller  than  that  of  Lag.33Si6026  when  doping  M2+  at 
the  La  site.  Furthermore,  samples  with  higher  dopant  levels  result 
in  lower  o  values.  For  instance,  the  o  values  of  La8M2Si6  026  (open 
circles)  are  lower  than  those  of  La8.67MSi6026  (shaded  circles).  Such 
trends  are  in  a  good  agreement  with  reported  observations  [3] 
which  revealed  that  the  magnitude  of  the  change  in  o  to  be  strongly 
dependent  on  the  size  of  the  dopant  cation,  and  that  small  dopants 
at  the  La  site  are  detrimental  to  a.  This  study  showed  that  ionic 
size  may  not  be  the  sole  determining  factor.  For  instance,  the  ionic 
radii  (A)  of  Mg2+,  Ca2+,  Sr2+  and  Ba2+  are  0.72,  1,  1.18  and  1.35, 
respectively,  while  the  ionic  radius  of  La3+  (=1.03)  is  in  between. 
The  present  study  reveals  that  another  factor,  the  electronegativ¬ 
ity,  also  plays  an  important  role.  In  fact,  the  coupling  effect  of  ionic 
radius  and  electronegativity  determines  the  resulting  ionic  con¬ 
ductivity  via  the  relative  Coulomb  energy.  The  calculation  shows 
that  the  Coulomb  energy  of  Lag.33Si6026  is  0.191,  while  the  resul¬ 
tant  relative  Coulomb  energies  of  La8M2Si6026  are  0.172,  0.177, 
0.180  and  0.183  for  M  =  Mg2+,  Ca2+,  Sr2+  and  Ba2+,  respectively,  that 
is  all  values  are  below  0.191.  Similarly,  the  Coulomb  energies  of 
La8.67MSi6026  (M  =  Mg2+,  Ca2+,  Sr2+  and  Ba2+)  are  also  less  than 
0.191. 

Similar  to  the  first  scenario,  two  levels  of  dopant  contents 
were  investigated  for  the  second  and  third  scenarios,  respectively. 
The  two  dopant  levels  are  represented  by  the  general  formula  of 
La8.33MSi6026  and  La7.33M2Si6  026-  The  calculated  a  values  of  the 
second  scenario  are  presented  in  Fig.  5b  (shaded  and  open  circles). 
The  behavior  of  this  group  of  trivalent  cations  is  generally  similar 
to  that  of  the  divalent  cations  in  the  first  scenario,  except  that  Ga3+ 
deviates  slightly  from  the  trend.  The  calculated  o  values  of  the  third 
scenario  are  shown  in  Fig.  5c  (shaded  and  open  circles).  As  seen 
in  Fig.  5b  and  c,  all  the  calculated  o  values  are  lower  than  that  of 
Lag.33Si6026,  and  the  larger  doping  content  leads  to  a  more  signifi¬ 
cant  decrease  in  o.  These  results  can  be  explained  from  the  relative 
Coulomb  energies  of  the  samples  in  the  two  scenarios.  For  instance, 
the  relative  Coulomb  energies  of  La8.33MSi6026  (M  =  Bi3+,  Nd3+,  Pr3+ 
and  Ce3+)  are  0.1862,  0.1904,  0.1905  and  0.1908,  respectively,  all 
of  which  are  smaller  than  that  of  Lag.33Si6  026-  A  further  examina¬ 
tion  of  the  third  scenario  shows  that  the  same  trends  as  those  of 
the  first  scenario  are  observed  in  the  samples  doped  with  Nd3+, 
Pr3+  and  Ce3+.  It  is  noted  that  the  ionic  radii  of  these  three  cations 
increase  (RNd  <RPr  <Rcel  while  their  electronegativities  decrease 
(ENd  >  Epr  >  ECe  ).  In  the  case  of  dopant  Bi3+,  its  ionic  radius  is  the  same 
as  La3+ ,  but  its  electronegativity  ( 1 .9)  is  significantly  larger  than  that 
of  La3+  (1.1),  which  leads  to  the  significant  lower  relative  Coulomb 
energy  of  La8.33BiSi6026  (0.1 862).  This  explains  the  dramatic  reduc¬ 
tion  in  a  of  Bi-doped  oxyapatites  reported  in  [3].  Overall,  doping  at 
the  La  site  leads  to  a  reduction  of  ionic  conductivity,  but  the  dopants 
with  large  ionic  radii  and  small  electronegativities  at  the  La  site 
generally  result  in  less  decrease  in  ionic  conductivity. 

3.2.2.  Doping  at  the  Si  site 

In  the  same  manner  as  Section  3.2.1,  the  effects  of  doping  at 
the  Si  site  of  Lag.33Si6026  were  studied  for  five  simulated  scenarios. 


In  particular,  Scenarios  (l)-(3)  focus  on  the  same  cations  as  those 
doped  on  the  La  site  (see  Section  3.2.1 ).  Similarly,  two  different  lev¬ 
els  of  dopant  contents  are  calculated  for  each  of  the  three  scenarios 
and  the  results  are  plotted  in  Fig.  5a-c,  respectively,  in  comparison 
with  the  results  of  Section  3.2.1 .  Doping  tetravalent  cations  at  the  Si 
site  is  studied  in  Scenario  (4),  and  doping  trivalent  and  tetravalent 
cations  standing  in  the  4th  and  5th  period  in  the  periodic  table  is 
studied  in  Scenario  (5). 

The  results  of  Scenario  (1 )  are  shown  in  Fig.  5a  with  two  general 
formulae  of  Lai0Si5lVIO26  (shaded  diamonds)  and  Lag.67Si5.5M0.5O26 
(open  diamonds).  Doping  at  the  Si  site  with  Group  IIA  cations  shows 
a  steady  boost  in  ionic  conductivity,  a ,  and  the  higher  level  of 
dopants  results  in  a  larger  increase  in  cr,  which  is  similar  to  the 
trends  of  doping  at  the  La  site.  In  contrast  to  doping  at  the  La 
site,  however,  all  the  calculated  o  values  are  larger  than  that  of 
Lag.33Si6026  for  all  levels  of  dopants  at  the  Si  site.  This  finding  also 
agrees  well  with  the  reported  observation  in  [3]. 

The  results  of  Scenario  (2)  are  shown  in  Fig.  5b  with 
two  general  formulae  of  Lag.67Si5M026  (shaded  diamonds)  and 
Lag.5Si5.5M0.5O26  (open  diamonds).  It  is  found  that  doping  Group 
IIIA  cations  at  the  Si  site  results  in  trends  that  are  opposite  to  those 
of  doping  at  the  La  site,  and  the  cr  values  of  both  levels  of  dopant 
contents  are  smaller  than  that  of  Lag  33Si6026-  On  the  other  hand, 
as  shown  in  Fig.  5c  (Scenario  3),  doping  rare  earths  of  Nd3+,  Pr3+ 
and  Ce3+  leads  to  a  boost  in  o.  The  substitution  of  Bi3+  at  the  Si  site, 
however,  lowers  the  ionic  conductivity  which  is  similar  to  that  at 
the  La  site. 

In  Scenario  (4),  the  substitution  of  isovalent  ions  at  the  Si  site 
with  the  general  formula  of  Lag.33Si5M026  is  examined  and  the 
results  are  presented  in  Fig.  6.  Two  opposite  effects  are  observed. 
Doping  W4+,  V4+,  Nb4+,  Ti4+,  Ta4+,  Zr4+  and  Uf4+  leads  to  an  increase 
in  cr,  where  as  doping  Se4+,  Rh4+,  Pt4+,  Pd4+,  Sn4+,  Ge4+  and  Pb4+ 
reduces  a.  While  the  ionic  radii  of  all  these  cations  in  the  two  groups 
are  nearly  equal  to  each  other,  the  electronegativities  of  the  cations 
of  the  former  group  are  obviously  smaller  than  those  of  the  latter 
group.  The  resultant  relative  Coulomb  energies  of  the  former  group 
are  therefore  larger  and  those  of  the  latter  group  are  smaller  than 
the  Coulomb  energy  of  Lag.33Si6  026- 

The  effects  of  doping  at  the  Si  site  with  trivalent  and  tetrava¬ 
lent  cations  are  compared  further  in  Scenario  (5).  Two  types  of 
oxyapatites  are  evaluated  with  general  formulae  of  Lag.67Si5M026 
(M  =  Ni3+,  Co3+,  Fe3+,  Pd3+,  Rh3+,  Ru3+)  and  Lag^SisMC^e  (M  =  V4+, 
Ti4+,  Nb4+,  Zr4+).  The  results  are  presented  in  Fig.  7.  Apparently,  the 
ionic  conductivities  of  oxyapatites  doped  with  tetravalent  cations 
are  noticeably  higher  than  those  with  trivalent  cations.  It  is  noted 
that  the  electronegativities  of  V,  Ti,  Nb  and  Zr  are  significantly 
smaller  than  those  of  Ni,  Co,  Fe,  Pd,  Rh  and  Ru,  while  the  ionic 
radii  of  the  tetravalent  cations  are  slightly  larger  or  equivalent 
to  those  of  the  trivalent  cations.  The  coupling  effect  of  these  two 
properties  leads  to  higher  ionic  conductivities  in  oxyapatites  doped 
with  tetravalent  elements.  In  summary,  the  calculation  shows  that 
doping  at  the  Si  site  with  cations  of  larger  ionic  radii  and  smaller 
electronegativities  generally  results  in  more  significant  increases 
in  ionic  conductivity. 

While  the  predicted  results  suggest  that  substitution  of  Ti  for  Si 
improves  the  ionic  conductivity,  it  seems  to  contradict  the  observa¬ 
tion  reported  in  [3  ].  The  reported  low  ionic  conductivity  of  Ti-doped 
oxyapatite  could  be  related  to  the  presence  of  secondary  phase  [28]. 

3.2.3.  Effect  of  ionic  radius  and  electronegativity 

It  has  been  shown  that  the  ionic  conductivity  of  an  oxyapatite 
is  determined  by  the  coupling  effect  of  the  ionic  radii  and  the  elec¬ 
tronegativities  of  the  constituents  via  the  relative  Coulomb  energy. 
To  have  a  better  insight  into  the  roles  of  the  two  types  of  fundamen¬ 
tal  property,  simulation  studies  are  performed  by  systematically 
varying  the  ionic  radii  of  dopants  in  the  range  of  0.2- 1.4  A  and 
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Ionic  Radius 

Fig.  8.  Contour  plots  of  ionic  conductivity  (Scm  1 )  variation  affected  by  ionic  radius  and  electronegativity,  (a)  Las.67MSi6026  (M  =  2+),  (b)  La8.33MSi6  026  (M  =  3+),  (c) 
La10Si5MO26  (M  =  2+),  (d)  Lag^SisMCfee  (M  =  3+),  (e)  Lag^SisMOze  (M  =  4+). 


electronegativities  in  0.7-4.  The  study  covered  dopants  of  valence 
2+  and  3+  at  the  La  and  Si  sites,  respectively,  and  valence  4+ 
at  the  Si  site.  By  substituting  one  cation  (denoted  as  M)  into  La 
or  Si,  five  types  of  pseudo  oxyapatite,  namely,  (a)  La867MSi6026 
(M  =  2+),  (b)  La8.33MSi6026  (M  =  3+),  (c)  La10Si5MO26  (M  =  2+), 


(d)  Lag  67Si5M025  (M  =  3+)  and  (e)  Lag33Si5M025  (M  =  4+)  are 
generated. 

Five  scenarios  that  are  associated  with  the  five  types  of  oxya- 
patites  are  obtained  and  the  contour  plots  are  given  in  Fig.  8.  A 
number  of  bands  are  observed  in  each  plot  and  correspond  to  dif- 
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Table  3 

Predicted  ionic  conductivity  a  for  new  oxyapatites. 


Type 

Oxyapatite 

cr  (Scm-1) 

A9.33B6O26  (A3+,  B4+) 

Ce9.33  ( Hf,Ta,Ti,Zr  )6  026 

1.1  x  10-3-8.7x  10-3 

La9.33(Hf,Nb,Ta,Ti,V,Zr)6026 

1.7  x  10_3-4.7  x  10-3 

Nd9.33(Hf,Zr)6026 

1.1  xl0-3-1.5xl0-3 

Pr9.33  (Hf,Zr)6026 

2.1  x  10“3-2.9x  10-3 

AioB15B2026  (A3+,  Bl4+,  B22+) 

CeioHf5(Ba,Ca,Sr)026 

1.7  x  10-2-2.7  x  10-2 

CeioZr5(Ba,Ca,Sr)026 

1.3x10-2-2.0x10-2 

Laio  Hfs  ( Ba,Ca,Sr  )026 

6.8x10-2-1.1  xlO-1 

LaioHf5(Co,Cu,Fe,Mg,Ni,V,Zr)026 

1.1  x  10-2-3.4  x  10-2 

LaioTa5(Ba,Ca,Sr)026 

1.2x10-2-2.0x10-2 

LaioTi5(Ba,Ca,Sr)026 

1.2x10-2-1.9x10-2 

Lai0Zr5(Ba,Ca,Sr)O26 

5.1  x  10-2-8.1  x  10-2 

La10Zr5(Fe,Mg,V,Zr)O26 

1.1  x  10-2-2.5x10-2 

A2BI4B22O26  (A3+,  Bl4+,  B23+) 

Ce2Hf4(Ce,Dy,Eu,Gd,Ho,La,Lu,Nd,Pr,Sm,Tb,Y,Yb)2026 

1.0  x  10-2-3.7  x  10-2 

Ce2Zr4(Ce,Eu,La,Lu,Nd,Pr,Sm,Tb,  Yb)2026 

1.1  x  10-2-2.9  x  10-2 

La2Hf4(Ce,Eu,La,Lu,Tb,Yb)2026 

8.2  x  10-2-1.5  xlO”1 

La2Hf4(Dy,Er,Gd,Ho,Nd,Pr,Sm,Y)2026 

3.9x10-2-7.6x10-2 

La2Nb4(Eu,La,Lu,Tb,Yb)2026 

1.1  x  10-2-1.9  x  10-2 

La2Si4(Eu,La,Lu)2026 

1.0  x  10-2-1.6  x  10-2 

La2Ta4(Ce,Dy,Eu,Gd,Ho,La,Lu,Nd,Pr,  Sm,Tb,  Y,Yb)2026 

1.0x10-2-3.9x10-2 

La2Ti4(Ce,Dy,Eu,Gd,Ho,La,Lu,Nd,Pr,  Sm,Tb,Y,Yb)2026 

1.1  x  10-2-3.9  x  10-2 

La2V4(Ce,Eu,La,Lu,Nd,Pr,Tb,Yb)2026 

1.0  x  10-2-2.1  x  10-2 

La2W4Lu2026 

1.0  x  10-2 

La2Zr4Lu2026 

1.2  x  10-1 

La2Zr4(Ce,Dy,Er,Eu,Gd,Ho,La,Nd,Pr,Sm,Tb,Y,Yb)2026 

3.0x10-2-7.3x10-2 

Pr2Hf4Lu2026 

1.2x10-2 

A9.33BI5B2O26  (A3+,  Bl4+,  B24+) 

La9.33Hf5(Ge,Nb,Pb,Si,Sn,Ta,Ti,V,W,  Zr)026 

1.1  x  10-2-3.2  x  10-2 

La9.33Zr5(Hf,Nb,Pb,Si,Ta,Ti,V,W)026 

1.0x10-2-2.6x10-2 

ferent  levels  of  ionic  conductivity.  The  higher  ionic  conductivities 
are  generally  located  at  the  lower  right  corners  of  the  plots  and 
are  associated  with  larger  radii  and  lower  electronegativities.  The 
ionic  conductivities  a  are  essentially  below  1  x  10-3  Son-1  when 
doping  divalent  cations  at  the  La  site  (Fig.  8a).  By  contrast,  doping 
trivalent  cations  leads  to  an  increase  in  o  if  the  ionic  radii  are  larger 
than  1 A  and  electronegativities  are  smaller  than  1.3  (Fig.  8b).  For 
doping  at  the  Si  site  (Fig.  8c-e),  all  the  dopants  of  valence  2+,  3+ 
and  4+  enhance  the  ionic  conductivities  if  their  electronegativities 
are  smaller  than  1.3.  In  this  case,  the  ionic  radius  is  not  a  dominant 
factor,  although  a  cation  with  larger  ionic  radius  results  in  a  higher 
cr. 

3.3.  Prediction  of  high  ionic  conductivity 

In  order  to  offer  a  guideline  for  the  search  of  new  oxyapatites 
with  high  ionic  conductivities,  the  ionic  conductivities  for  a  wide 
range  of  pseudo  oxyapatites  are  predicted  by  using  the  developed 
model.  The  new  oxyapatites  are  virtually  generated  by  numerically 
substituting  a  wide  range  of  cations  with  valence  of  2+,  3+,  and  4+ 
into  the  La  and  Si  sites,  respectively.  As  mentioned  above,  doping 
divalent  cations  at  the  La  site  is  detrimental,  while  doping  trivalent 
cations  can  increase  the  ionic  conductivity  only  if  the  ionic  radius 
is  larger  than  1 A  and  the  electronegativity  is  smaller  than  1.3  (as 
indicated  in  Fig.  8).  Unfortunately,  only  a  few  trivalent  cations  meet 
such  requirements.  On  the  other  hand,  doping  at  the  Si  site  is  gen¬ 
erally  improves  the  conductivity.  Therefore,  four  types  of  pseudo 
oxygen  stoichiometric  oxyapatite  are  generated  with  the  general 
formula  as  (1)  A9.33B6026  (A3+,  B4+),  (2)  Ai0B15B2O26  (A3+,  Bl4+, 
B22+),  (3)  A10Bl4B22O26  (A3+,  Bl4+,  B23+),  and  (4)  A9.33B15B2026 
(A3+,  Bl4+,  B24+).  In  particular,  the  following  cations  are  covered. 

•  A-site:  Ce3+,  Dy3+,  Er3+,  Eu3+,  Gd3+,  Ho3+,  La3+,  Lu3+,  Nd3+,  Pr3+, 
Sm3+,  Tb3+,  Y3+,  Yb3+,  B3+,  Al3+,  Ga3+,  In3+,  Tl3+,  Bi3+. 

•  Bl-site:  Ge4+,  Se4+,  Si4+,  Sn4+,Ta4+,Ti4+,  V4+,  W4+,  Zr4+,  Nb4+,  Hf4+, 
Rh4+,  Pt4+,  Pd4+,  Pb4+. 


•  B2-site:  (a)  divalent:  Ba2+,  Ca2+,  Co2+,  Cr2+,  Cu2+,  Fe2+, 
Mg2+,  Mn2+,  Mo2+,  Ni2+,  Pd2+,  Sr2+,  V2+,  Zn2+;  (b)  trivalent: 
same  cations  as  A-site;  and  (c)  tetravalent:  same  cations  as 
Bl-site. 

Since  the  number  of  such  pseudo  oxyapatites  is  very  large, 
only  those  with  calculated  high  ionic  conductivities  (at  500  °C)  are 
selected  and  listed  in  Table  3.  The  ionic  conductivities  of  type-1 
(A9i33B6  026)  are  the  lowest  (below  1  x  10-2  S  cm-1)  among  the  four 
types  of  oxyapatite.  By  contrast,  the  ionic  conductivities  of  the  other 
types  are  generally  of  the  order  of  10-2  S  cm-1,  and  some  of  them 
are  as  high  as  10_1  Son-1.  It  is  also  found  that  La  or  Ce  should  be 
predominate  at  the  A  site  in  order  to  obtain  high  ionic  conductiv¬ 
ities.  On  the  other  hand,  the  substitution  of  Zr  or  Hf  at  the  B  site 
leads  to  the  highest  a.  To  investigate  the  effect  of  doping  levels  of 
Zr  at  the  Si  site,  the  ionic  conductivity  of  Lag^Sig-yZryO^  is  calcu¬ 
lated  withy  varying  in  the  range  of  0-6.  The  results  are  presented  in 
Fig.  9.  The  ionic  conductivity  increases  steadily  with  increasing  the 
Zr4+  level.  Furthermore,  a  significant  increase  in  ionic  conductivity 
occurs  only  when  y  >  4. 


y 

Fig.  9.  Predicted  ionic  conductivity  of  La9.33Si6-yZry026  with  variation  in  doped  Zr 
content. 
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4.  Conclusions 

It  is  demonstrated  that  two  types  of  fundamental  chemical  prop¬ 
erty,  i.e.,  electro-negativities  and  ionic  radii  of  the  constituents, 
control  oxygen  ionic  conduction  in  oxyapatites.  The  two  proper¬ 
ties  have  been  used  to  represent  the  relative  charge  densities  and 
the  distances  between  the  charged  units,  respectively,  and  then  to 
formulate  the  relative  Coulomb  energy.  It  is  found  that  this  relative 
Coulomb  energy  is  linearly  correlated  to  the  oxygen  ionic  conduc¬ 
tivity  (in  logarithmic  form)  in  the  oxyapatite  systems.  By  applying 
the  concept  of  the  relative  Coulomb  energy,  most  experimentally 
observed  trends  can  be  well  rationalized. 

A  wide  range  of  new  oxygen  stoichiometric  oxyapatites 
have  been  predicted  with  ionic  conductivities  in  the  order  of 
10-2-10-1  Scm-1.  Among  them,  lanthanum-based  oxyapatities 
with  Zr4+  or  Hf4+  substituting  for  Si4+  have  the  highest  ionic  con¬ 
ductivities,  which  implies  that  this  family  of  oxyapatites  may  be 
promising  candidates  for  the  electrolyte  in  IT-SOFCs.  The  predic¬ 
tion  also  shows  that  the  substitution  of  Ta4+,  Ti4+,  Nb4+,  W4+  and 
V4+  for  Si4+  and  the  substitution  of  Ce3+  for  La3+  results  in  rela¬ 
tively  higher  ionic  conductivities.  Finally,  this  study  reveals  that 
doping  a  cation  with  large  ionic  radius  and  low  electro-negativity 
tends  to  increase  the  ionic  conductivity  of  the  oxyapatite.  While 
these  results  are  required  to  be  verified  by  experiments,  they  offer 
a  guideline  in  searching  for  promising  oxyapatites  with  high  ionic 
conductivities. 
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